All-optical Imprinting of Geometric Phases onto Matter Waves 
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Traditional optical phase imprinting of matter waves is of a dynamical nature. In this paper we 
show that both Abelian and non-Abelian geometric phases can be optically imprinted onto matter 
waves, yielding a number of interesting phenomena such as wavepacket re-directing and wavepacket 
splitting. In addition to their fundamental interest, our results open up new opportunities for robust 
optical control of matter waves. 
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When the Hamiltonian of a quantum system adiabati- 
cally changes along a path in its parameter space, a non- 
degenerate (degenerate) eigenstate will acquire, on top of 
an obvious dynamical phase, an Abelian (non-Abelian) 
geometric phase Q . The discovery of geometric phases 
has motivated many proposals on built-in fault-tolerant 
quantum logical operations [3]. Geometric phases also 
provide an elegant framework to understand the coupling 
between translational motion and internal degrees of free- 
dom. In particular, a slowly moving system subject to 
an external field naturally causes adiabatic changes in 
its internal Hamiltonian, giving rise to an Abelian or 
non-Abelian vector potential for the translational motion 
0, H, 0] • Interestingly, one can also engineer laser-matter 
interaction to synthesize effective vector potentials from 
geometric phases [7, 8, 9, 10]. 

One implicit assumption in geometric phase-based 
quantum gate studies is that the qubit systems are frozen 
in space. By contrast, here we examine the geometric 
phases acquired by a system coherently delocalized in 
space. In particular, we propose to optically imprint ge- 
ometric phases at each local space point onto delocalized 
matter waves of ultracold systems. We show that the 
imprinted geometric phases can alter the ensuing matter- 
wave propagation dramatically, thus establishing another 
interesting and potentially powerful means of optical con- 
trol. Because here the geometric phases are not induced 
by the translational motion itself, but by active manip- 
ulation of laser-matter interaction, the imprinting of ge- 
ometric phases onto matter waves results in a new type 
of coupling between internal and translational motions, 
which can be implemented experimentally. 

In generating matter-wave vortices with time- 
dependent magnetic fields, some early work did take ad- 
vantage of a spatially-dependent Abelian Berry phase 
[ll|. Our work however represents a major extension 
from the vortex context, by demonstrating the imprint- 
ing of both Abelian and non-Abelian geometric phases 
with all- optical methods. We shall show that imprinting 



geometric phases may re-direct atom wavepackets, split 
a wavepacket, or create interesting matter wave patterns. 
The all-optical geometric phase imprinting is insensitive 
to the dynamical details. Hence it is superior to con- 
ventional optical phase imprinting approach to matter 
wave engineering [12*, where the imprinted phase is 
proportional to the duration of optical imprinting and to 
the magnitude of an optical potential. 

Consider first a two-level atom in a plane-wave laser 
field. In the rotating wave approximation (RWA), the in- 
ternal Hamiltonian is given by i^RWA,2 = (^/2)(|1)(1| — 
|2)(2|) + (l]o/2)(e-^^^-^|l)(2|+e^k'-i^|2)(l|), where |1) and 
1 2) are the two internal states of the atom, is assumed 
to be real, is the wavevector of the laser field, A is the 
detuning from one-photon resonance, and R the coordi- 
nate of the atom. In the (|1), |2)) representation i^RWA,2 
has the eigenstate 
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X(A) 
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and an analogous eigenstate ) orthogonal to |(/>^), 

where x(A) = yj^l^ iV^o + + ^ is the normal- 
ization factor. Note that these internal eigenstates vary 
with the coordinate R. Within the adiabatic approxima- 
tion and in the absence of a trapping potential, the over- 
all eigenstate (internal plus translational) for the system, 
denoted |^^), is given by |^^) = |(/)^)e*^ "^, where k is 
a wavevector associated with the translational motion of 
the atom. 

We begin with a test case that is well-known in laser 
chirping control [14] (which realizes a Landau-Zener pro- 
cess). Here the detuning parameter A is adiabatically 
varied from Ai > to A2 < 0, with \Ai/Qo\ > 1 and 
|A2/l^o| ^ 1- Accordingly k^ changes from k^^ to k2, but 
with a fixed direction. The initial state is assumed to be 
|^+(A = Ai)) ^ (l,0)^e^^-^. Assuming that the trans- 
lational motion during the chirping process is negligible, 
which can be a good approximation for ultracold atoms. 
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one finds the following geometric phase for an arbitrary 
R, 
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where kj, denotes the laser wavevector when A = [15| , 
and throughout k represents a unit vector along the vec- 
tor k. Given that the evolving state remains in the 
system eigenstate, the dynamical phase is evidently R- 
independent. Imprinting the geometric phase /3(R) onto 
the adiabatic state |^+) in Eq. ([Tj) and neglecting the 
overall dynamical phase, one obtains the final state 



l*l(A2)> 



gi[/3(R)+k-R] 
X(A2) 



i(k^+k)-R 
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Equation (|3]) shows that the population between the two 
levels is adiabatically inverted and the atom momentum 
increases by Tiki. This result is expected, because (i) 
the chirping scenario is long known to induce the ab- 
sorption of one on-resonance photon and (ii) momentum 
conservation requires the atom momentum to increase by 
Tiki. But quite noteworthy, here the obvious consequence 
of momentum conservation is manifested as an effect of 
the spatially-dependent geometric phase /3(R). Note also 
that if we further change A along a reversed path until it 
reaches its initial value, then Eq. ([2]) gives a zero Berry 
phase and hence zero momentum change. This is again 
correct because the atom releases a photon during the 
reversed process. One can also carry out similar analy- 
sis if the chirped laser field possesses a nonzero orbital 
angular momentum. That is, the transfer of orbital an- 
gular momentum from photon to the matter wave (hence 
vortex generation) can also be interpreted as geometric 
phase imprinting. 

As a second Abelian example, we examine the sys- 
tem of i^RWA,2 for A = and with k^ being rotated 
adiabatically. We assume below both kj, and k are 
in the x-z plane. Let 7 be the angle by which kj, 
is rotated clockwise with respect to the x-axis. Then 
k^(7) = kl [cos(7)ea, - sin{j)ez] , where = |kj,|, Cx and 
Cz are the unit vectors along the x-axis and z-axis. For 
convenience we also define the polar coordinates (i?, 6) 
for R confined on the x-z plane, i.e., x = Rcos{0)^ and 
z = Rsm{0). Then, the eigenstate in Eq. ([1]) for A = 
can be rewritten as 
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,i/cJ.i?cos(6'+7) 



(4) 



Consider an initial state |^+) = |^+(7 = 0))e*^'^. As 
the adiabatic parameter 7 varies from to a, the acquired 
geometric phase at an arbitrary local point (i?, 0) is found 
to be 



[Rcos{0) - Rcos{0^a)]. (5) 



Neglecting the translational motion during the adia- 
batic process, the final state can be obtained 
by imprinting P{R,6) onto the adiabatic state |^^(7 = 
^j^gzk-R^ yielding 



|f+(a)) 



i[/3(i?,6')+k-R] 
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I / i/ct.i?sin(0+f)sin(f) \ 
I gifej,i?cos(0+-f ) cos(-f ) I ^ 
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Unlike in the above test case that apparently involves the 
absorption of one on-resonance photon, here neither the 
population on states |1) and |2) nor the atom mechani- 
cal momentum changes. However, the geometric phases 
P{R^6) imprinted onto the matter wave will induce dif- 
ferences in time evolution if we suddenly switch off the 
laser field, after which the matter waves associated with 
each of the two components of (a)) will evolve inde- 
pendently and freely. The momentum in free space for 
each component of (a)) can be found by taking the 
spatial derivative of that component. Without loss of 
generality we set k = 0. One then finds that upon a sud- 
den laser field switch-off, the state (a)) will split into 
two components: atoms with the internal state |2) will 
move in the direction of cos(a/2)ea; — sm{a/2)ez with a 
momentum cos(a/2), but the atoms on the state |1) 
will move in the direction of sm{a/2)ex^cos{a/2)ez with 
a momentum hkl sm{a / 2) . These motion characteristics 
are clearly a-dependent and intriguing. In particular, if 
there is no geometric phase imprinting, a = 0, the |1) 
component is static and the |2) component moves along 
Cx] but if a = TT, then the |2) component is static and 
the |1) component moves along Cx. This could be useful 
for filtering the internal states. 

We have also carried out numerical wavepacket simula- 
tions based on the Hamiltonian — + ^rwa,2 be- 
ing the mass of the atom) that undergoes slow changes in 
its parameter 7. In our simulations we replace the plane- 
wave factor exp(zk-R) of |^+) by a Gaussian wavepacket. 
Numerical results do confirm our predictions. 

We now discuss the imprintingof non- Abelian geo- 
metric phases. As in Refs. [iq|, llfil, we adopt the 
so-called "tripod scheme" of four-level atoms interact- 
ing with three on-resonance laser fields. The internal 
state is denoted as |n), n = — 3. Each of the three 
transitions |0) ^ |0) ^ |2), and |0) ^ |3) is cou- 
pled by one laser field. This coupling scheme can be 
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realized if states |1), |2), and |3) are degenerate mag- 
netic sub-levels and the three coupling fields have dif- 
ferent polarizations. Alternatively, states |1), |2), |3) 
are non-degenerate and the three coupling fields then 
carry different frequencies (same polarization is then al- 
lowed). For convenience we adopt the same configura- 
tion as in Ref. [13], where two laser beams are counter- 
propagating along the x-axis and the third laser beam 
is along the 2:-axis. The associated internal Hamiltonian 
under RWA is given by i^RWA,4 = Yln=i ^n|0)(n| -\- h.c.^ 
with Qi = Qo^HO/V^e-'^'r-^, ^2 = l^osin(0/y2e'^''^, 
and Qs = cos{^)e'^^^^ , where the parameter ^ is set 
to satisfy cos(^) = \/2 — 1. Note that even if the three 
lasers carry different frequencies, their effective kl in the 
x-z plane can still be the same by titling the laser beams 
out of the x-z plane. 

The Hamiltonian i^RWA,4 has two degenerate and 
spatially-dependent dark (null-eigenvalue) states |l^i(2)), 
which are given by [17], 



\Dr) = (|i)-|2))e- 
11^2) = [cos(0 + |2)) /V2 - sin(0|3)] e^' ^(7) 

where f^' = kl[l - cos(0], |i) = |l)e*^' (^+^), and |2) = 
|2^e-*^r(3^-^). Clearly, within the dark-state subspace the 
dynamical phase is always zero. So we focus on the time 
evolution in this dark subspace. Any state therein can be 
expanded as ci\Di) + C2|I^2)- This expansion hence de- 
fines a dark state representation. In this representation, 
the mechanical momentum operator becomes pjl, ll8| 



(8) 



where (Jx^z are Pauli matrices, n = cos(^)/c[, and V repre- 
sents the gradient in the dark-state representation. As in 
Ref. [17], we also introduce an additional constant shift 
to state 1 3), denoted by a matrix operator F^, such that 



the total effective Hamiltonian becomes Hf^ — ^ 
The eigenstates of are 



1^ 



1 T ^e*^' 
-i ± e^^^ 



^ik R 



(9) 



where is the angle between the x-axis and the atom 
wavevector k in Eq. ([9|). According to Eq. (|8|), it is 
straightforward to show that the eigenstates |^^'^) have 
the momentum (k ± A>:k)fi. 

Consider then two simple scenarios of laser field ma- 
nipulation. In the first scenario, we adiabatically move 
all the laser beams along the negative z axis. The non- 
Abelian geometric phases thus induced are determined 
by 
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FIG. 1: (Color Online) Numerical simulation of wavepacket 
splitting induced by imprinting non-Abelian geometric 
phases. Here t is in units of m/{hHi^); x and z are in units of 
1/k. See the text for details. 



Because the matrix on the right side of Eq. (p!Q|) turns 
out to be independent of R, one may naively conclude 
that the non-Abelian geometric phases induced here will 
not affect matter wave propagation. But this is incor- 
rect. To see this let us assume the initial state to be the 
1^^'") state in Eq. ^ with cpk = 0. This initial state 
has a mechanical momentum h{k — n) along the x-axis, 
where k = |k|. Let the final state, a function of the dis- 
placement dz along — be {dz))- It is enlightening to 
consider a specific case, e.g., dz = f;^- Neglecting matter 
wave propagation during the period of geometric phase 
imprinting, Eq. (p!Q|) then gives 



1 + i f 1 
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Akx 



.(11) 



The two components on the right side of Eq. (pT]) rep- 
resent a superposition of the two eigenstates |^^'+) and 
|^^'~) in Eq. ([9|). Because these two eigenstates pos- 
sess different mechanical momenta {k^hz)h and {k — K)h, 
we thus have the result that the non-Abelian geometric 
phases generated here can split matter waves. 

This result is also verified by our numerical simula- 
tions based on the full Hamiltonian +-^rwa,4 + '^- 
In particular, in our simulations the initial state is cho- 
sen as a Gaussian wavepacket instead of a plane wave 
considered in Eq. (|9]). As an example we choose 
Qo/h = 1200hf<i^ /m. The duration of the imprinting 
process is chosen to be 0.063m/ {Hk'^), (In real units, for 
m = 10~^^ Kg, ~ 10^ m~^, the duration is ~ 60 /iS and 
Qo/h ^ 10^ Hz. The cold-atom version of the Zitterbe- 
wegung oscillation [19] is negligible for these parameters). 
Figure 1 depicts the simulation results for (pj. = 0, k = hz. 
At t = we have initiated the first non-Abelian geomet- 
ric phase imprinting for d^ = 7r/(4/^). As seen in Figs. 
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1(b)- (c), an initial wavepacket of zero group velocity in- 
deed splits into two parts: one part holds a zero group 
velocity and the other part moves along the x axis at an 
average mechanical momentum of 2?in. Then, at t = 20, 
we move all the lasers up by dz = 7r/(4A>:) to induce a 
second imprinting process. Figure 3(d) shows that this 
further splits the wavepacket into four copies of equal 
weights, a result also derivable from Eq. (p!Q|) . Choosing 
a different at t = 20 for the second imprinting process, 
we may even exchange the momenta between the two sep- 
arated sub-wavepackets. Extending the strategy here, in 
principle one can split a wavepacket into 2^ copies with 
an arbitrary integer n. This should be useful for atom 
optics applications such as atom interferometry. Other 
simulation results (not shown) also indicate that so long 
as VLq is sufficiently large, then on one hand the matter 
wave propagation during the adiabatic processes can be 
made negligible and on the other hand the adiabaticity 
in the internal motion is well maintained. 

In the second scenario, we slowly rotate all the three 
laser beams clockwise in the x-z plane. The non-Abelian 
geomertic phase at the local point {R^O) is determined 

by 

. / ci \ _ / i{DM^\D^) ^{D^\^\D,) \ ( c^\ 
'<^^\C2) \^{D^\^\D^) i{D^\^\D^) ) \ c2 J 

y - sm[0 + 7) - cos[0 + 7) y \ C2 y 

where < 7 < a is the angle of rotation that serves 
as the adiabatic parameter. An analytical solution to 
Eq. (p!2|) is not found and hence we solve it computation- 
ally. After obtaining the numerical non-Abelian geomet- 
ric phases for each point (i?, 0) from Eq. (p!2]) . we imprint 
the spatially dependent geometric phases onto the system 

fF 

at t = and then let it evolve under = (also with 
the rotated laser fields). For the example shown in Fig. 2, 
the initial state is the state |^^'~) in Eq. (|9]) multiplied 
by a Gaussian profile, with k = n and (fk = — 7r/2. This 
initial state hence has a zero group velocity according to 
Eqs. dH]) and (|9|). As shown in Fig. 2 for a = 37r/2, 
the phase imprinting substantially impacts the ensuing 
evolution. The evolving wavepacket first displays inter- 
esting interference patterns, followed by a certain degree 
of self-focusing, and later it starts to display some moon- 
like patterns. As expected, the patterns in real space at 
later times begin to resemble the t = wavepacket pat- 
tern plotted as a function of and Pj^, namely, the 
mechanical momenta along the x-axis and z-axis. Re- 
markably, even if we let a = 27r, i.e., with all laser fields 
returning to their original configuration, similar results 
can still be obtained. We have checked that essentially 
the same results can be obtained using the full Hamil- 
tonian — — h i^RWA,4 + Vs. However, we find that 
in order to fully neglect matter wave propagation during 
the imprinting, the duration of the adiabatic process here 
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FIG. 2: (Color Online) Simulation of the wavepacket dynam- 
ics after imprinting non-Abelian geometric phases via adia- 
batic rotation of the three laser beams in a tripod scheme. 
The wavepacket plotted in momentum space (top panel, in 
units of lity) will remain the same for any t > 0. Length units 
(for the bottom four panels in x-z space) are the same as in 
Fig. 1. 

should be about two orders of magnitude shorter than in 
the first non-Abelian example. Hence VLq must be about 
two orders of larger. 

In summary, we have demonstrated the concept of all- 
optical imprinting of both Abelian and non-Abelian ge- 
ometric phases onto matter waves. The adiabatic paths 
considered here are among the simplest. Exploring more 
complicated adiabatic paths might offer extensive laser 
control of matter wave propagation. Extending this work 
to the imprinting of non-adiabatic geometric phases is 
also of considerable interest. 

This work was supported by WBS grant No. R-710- 
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Quantum Computation", and by the NUS "YIA" fund 
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